Lattice QCD at finite isospin density. 

J. B. Kogut a , D. K. Sinclair 15 * 

a Physics Department, University of Illinois, 1110 West Green Street, Urbana, IL 61801, USA 
b HEP Division, Argonne National Laboratory, 9700 South Cass Avenue, Argonne, IL 60439, USA 

We have simulated QCD at a finite chemical potential fii for isospin (J3) to probe part of the phase diagram 
for nuclear matter. Preliminary results suggest that for fii > p c , this theory forms a charged pion condensate 
which spontaneously breaks 73, and the isospin density is non zero. 



1. Introduction 

Nuclear matter has finite baryon number den- 
sity, finite (negative) isospin density and, under 
certain conditions, finite (negative) strangeness 
density. Since having a finite chemical potential 
for baryon number makes the fermion determi- 
nant complex and simulations intractable, we re- 
strict ourselves to zero baryon-number chemical 
potential. QCD at finite chemical potential \ii 
for isospin (J3) and zero chemical potential for 
baryon number, has a real positive fermion de- 
terminant and can thus be simulated. This de- 
scribes a surface in the phase diagram for nuclear 
matter. This theory turns out to be very similar 
to 2-colour QCD at finite quark-number density 
which we are also studying |]||. We will also 
use various devices to include the effects of finite 
chemical potential /is for strangeness. 

QCD at finite chemical potential fii for isospin 
has been studied by Son and Stephanov Q| us- 
ing effective (chiral) lagrangians. This analy- 
sis indicates is that it undergoes a transition at 
= /i c = m w to a state with a charged pion 
condensate, while the orthogonal charged pion 
becomes a true Goldstone boson. This conden- 
sate breaks I3 and parity spontaneously. 

Section 2 describes our simulations of 2 flavour 
QCD with a finite /ij. A brief discussion of how 
to include a finite us is included. Section 3 gives 
our conclusions. 



*Talk presented by D. K. Sinclair at Lattice2001, Berlin, 
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2. Lattice QCD at finite isospin chemical 
potential 

The lattice fermion action (staggered quarks) 
for QCD with isospin chemical potential, /i/ is 



(1) 



where is the standard staggered _pwith links 
in the +t direction multiplied by and those 
in the —t direction multiplied by e~2^. The 
term proportional to A, which explicitly breaks 
the third component of isospin, I3, is introduced 
to allow us to observe spontaneous breaking of 
this symmetry from finite lattice simulations. We 
will be interested in the limit A — > 0. The Dirac 
operator is 
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whose determinant 



detM = detL4U + A 2 ] 
where 

A = Ip(iii) + m. 
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Since this determinant is strictly positive, stan- 
dard simulation methods work. We use hybrid 
molecular dynamics simulations to tune the num- 
ber of flavours down from 8 to 2. 
We measure the chiral condensate 



(ipip) (xx), 
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the isospin condensate 
and the isospin density 
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At /ij = A = 0, there are 4 pseudo-Goldstone 
pions corresponding to 



(8) 
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Only 2 of these (those with i = 1,2) correspond 
to pseudo-Goldstone pions of the Nf =2 theory 
- the third Nf — 2 pion is a massive 3-link ob- 
ject at finite lattice spacing. However, the other 
2 Goldstone bosons should mimic the behaviour 
of the third Nf = 2 pion (710) in the continuum 
limit and are thus useful to consider. As fii in- 
creases one expects the i = 1 mass to fall, while 
the i = 2 mass should rise. The other 2 masses 
are expected to remain unchanged for small 
One expects a phase transition at fii = /i c ~ m n , 
above which isospin^) is spontaneously broken 
and the i = 1 state becomes a true Goldstone bo- 
son. This is only part of the story. Once j'g gains 
an expectation value, scalar mesons are also can- 
didate pseudo-Goldstone bosons. 

We are simulating this theory at (3 = 5.2, m = 
0.025,0.05, A = 0.005,0.0025 on an 8 4 lattice. 
We are also performing simulations at m = 0.05, 
A = 0.005 on an 8 3 x 4 lattice to study the finite 
temperature transition. 

The pion condensate is shown in figure |This 
strongly suggests that there is some fii = fi c 
above which isospin(/ 3 ) is spontaneously broken 
by a charged pion condensate. However, we will 
need to finish our simulations at A = 0.0025 and 
extrapolate to A = to validate this observation. 

The chiral condensate shown in figure |^ is ap- 
proximately constant for \ii < fj, c . For \ii > fj, c 
it starts to fall, indicating that the condensate is 
rotating from the chiral to the isospin-breaking 
direction — it does not appear to be a simple 
rotation which is suggested by the tree-level ef- 
fective Lagrangian approach, however. 
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Figure 1. The pion condensate as a function of 



The isospin density, shown in figure || is con- 
sistent with zero for /ij < fi c and commences to 
rise for fij > fi c . It reaches its saturation value 
(3) for m - 2. 

We will measure the spectrum of Goldstone and 
pseudo-Goldstone mesons on a larger (12 3 x 24) 
lattice. In addition we will store configurations 
for offline calculations of other hadron spectra 
(N, A, p...) which are also expected to exhibit 
interesting behaviour. 

3. Conclusions 

QCD at finite chemical potential for isospin(73) 
appears to spontaneously break isospin with a 
charged pion condensate, for \ij large enough. We 
need another A to verify this and to show that 
He > 0. We will repeat these simulations on a 
larger (12 3 x 24) lattice to measure the spectrum, 
in particular the Goldstone boson. Measurement 
of the instanton distribution is of interest. Is 
there any additional structure at high fij, where 
the effective Lagrangian approach is suspect, per- 
haps related to baryon thresholds? What we re- 
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Figure 2. The chiral condensate as a function of 
Hi. 



Figure 3. The density of the 3rd component of 
isospin as a function of /zj. 



ally would like to know is do such pion conden- 
sates form at finite baryon-number density (i.e. 
in nuclear matter) for \i\ large enough? 

Preliminary indications from our finite temper- 
ature (8 3 x 4) simulations are that the evaporation 
of the pion condensate for large enough tempera- 
ture occurs at a first order transition. This needs 
further study. 

Simulations with finite chemical potentials for 
isospin and strangeness will study the competi- 
tion between pion and kaon condensation [Q Two 
methods will be used to avoid a complex deter- 
minant. The first is to use a partially quenched 
approach where only the u and d quarks are dy- 
namical. The second is to include a c quark 
with m c = m s and a chemical potential for the 
3rd component of the SU(2) flavour group which 
mixes c and s. 

QCD at finite isospin, and QCD at finite 
isospin and strangeness should give sensible 
quenched results. In fact our earlier work on 
quenched QCD at finite baryon number density 
H can be reinterpreted as quenched QCD at fi- 
nite isospin density. 
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